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We would like to pursue a mathematical approach to turbulence that is 
able to predict bulk average quantities of experimental and practical signif- 
icance. Among such bulk average quantities, the Nusselt number N, the 
enhanced bulk average heat transfer due to convection, is perhaps the sim- 
plest physical objective of rigorous mathematical study of turbulence. 

The question we wish to address here is: what are the optimal bounds 
on N as a function of Rayleigh number Rl What is the effect of the Prandtl 
number? Theoretical results ([Q, f2|, [[|) predict a maximal ultimate behav- 
ior of N ~ v^R. The experimental findings indicate however that 

iV~ R q 

where the reported values for q belong approximately to the interval [|, |] 
for large R. The exponents 2/7 and 1/3 have been discussed by several 
authors, ([Q, ||). The recent results of (0) favor an exponent of 3/10 with 
logarithmic corrections. We will describe below some very simple rigorous 
results. 

The mathematical formulation of the problem is based on the three di- 
mensional Boussinesq equations for Rayleigh-Benard convection (||). These 
are a system of equations coupling the three dimensional Navier-Stokes equa- 
tions 

-t^ + u- Vu + Vp = crAu + crRTe z , V-u = (1) 
ut 
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to a heat advection-diffusion equation 



dT 

— + u ■ VT = AT. (2) 

at 

The five unknowns, incompressible velocity, u = (u,v,w), pressure p, and 
temperature T are functions of position x = (x, y, z) and time t. R is the 
Rayleigh number and a is the Prandtl number. e z = (0, 0, 1) is the unit vector 
in the vertical direction. For simplicity of exposition the equations have been 
non-dimensionalized: the vertical variable z is scaled so that it belongs to the 
interval [0, 1], the horizontal independent variables (x,y) belong to a square 
Q C R 2 of side L. The boundary conditions follows: all functions 

((u,v,w), p, T) are periodic in x and y with period L; u, v, and w vanish 
for z — 0,1, and the temperature obeys T — at z — 1, T — 1 at z — 0. We 
write 

|2 l I / I f(„ „ „\|2 
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\f(x,y,z)\ dxdydz 



for functions and vectors /. We use < • • • > for long time average: 

1 r* 



1 r 

(f) = lim sup - / f(s)ds. 
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We write / the horizontal average 

- 1 f L f L 

f = j2 J o J o f{x,y)dxdy. 

When a function depends on additional variables we write only the remaining 
variables after integration, so for instance 

|| Vu(-, t)\\ 2 = -^£ £ jT 1 | Vu(x, y, z,t)\ 2 dx dy dz. 

The Nusselt number is given in terms of the long time average of the 
vertical heat flux: 

N = l + (f b(z)dz^ (3) 

with 

1 fL fL 

b(z,t) = — w(x,y,z,t)T(x,y,z,t)dxdy = wT(z,t). (4) 
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A consequence of the equations of motion are two additional formulas for 
the Nusselt number: 

(||VT|| 2 ) = iV (5) 

and 

|Vu|| 2 ) = R(N-1). (6) 



The classical result of Howard, conditioned on on assumptions about 
statistical averages, is that N is bounded at very large Rayleigh numbers 
by a multiple of . The same kind of bound can be derived without any 
conditions (@). This bound is valid for all solutions, aspect ratios L and 
Prandtl numbers a and is also valid in a rotating frame, at arbitrary rotation 
speed. 

The system is not isotropic: the direction of gravity is singled out. Con- 
sider a function t(z) that satisfies r(0) = 1, r(l) = 0, and express the 
temperature as 

T(x,y,z,t) = T(z)+9(x,y,z,t). (7) 

The role of r is that of a convenient background (]nj, [0 ) that carries 
the inhomogeneous boundary conditions; thus 9 obeys the same homogeneous 
boundary conditions as the velocity. Note that, because r does not depend 
on x, y one has 

T{x, y, z, t) - T{z, t) = 9{x, y, z, t) - 9(z, t). (8) 

The equation obeyed by 9 is 

(d t + u ■ V - A) 9 = t" - wt' (9) 

where we used r' — The horizontal average of the vertical velocity van- 
ishes identically because of incompressibility 

w(z,t) = 0. 

Therefore the quantity b(z, t) can be written as 

b ( z ^) = [ J w ( x iyi z it) (6{ x ,V, z ,t) -0(z,t)) dxdy = w{T-T)(z,t) 



(10) 
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Multiplying the equation (|9|) by 9 and integrating we obtain 



n + (\\ve 



1 r'(z)b(z)dz\ + J\t\z)) 2 dz. 



We will choose the background profile r for simplicity to be a smooth profile 
concentrated in a boundary layer of width 5, 

' z s 



T[Z 



P 



with P(0) = 1 and P(s) = for s > 1. Using only elementary facts (fun- 
damental theorem of calculus, the boundary conditions and the Schwartz 
inequality) it is easy to see from flIU| ) that 

\b(z,t)\<z\\Vn(.,m-\MT-T)(-,t)\\ 

holds for any z. Let us define 

n = (\\V(T-T)f 

Note that, from the definition of n and @ it follows that 

n < N. 

From (0), (0) and (|) we obtain 

C 



with 



and 
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N < — + 2D5JR(N - l)vn 



'dP( s y 

ds 

dP(s) 



ds 



ds 



ds. 



Optimizing in 5 < 1 we get 

5- l = 



'^-{R(N-I)n}*+1 



and letting P —* s we obtain 

N < 2(Rn)^(N - 1)3 + 1 

Thus we have proved 



(12) 
(13) 



(14) 
(15) 



(16) 
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Theorem 1 Let 

n = (||V(T-T)|| 2X 

Then the Nusselt number [| ; \Q) for three dimensional Rayleigh-Benard 
convection satisfies 

N < 23 (Rn)s + 1. 

If one has no additional information then, using n < N in the inequality 
above we obtain the square-root bound 

N < + 1 

(The prefactor is not optimal. The search for optimal prefactors is better 
motivated for other systems, where the power law obtained rigorously co- 
incides with the one observed in experiments. When that is the case then 
the rigorous results can match experiments with remarkable accuracy ( [|T3|j ) ) . 
The exponent 1/3 (or anything less than 1/2 for that matter) have not been 
proven rigorously for the general system. The theorem above brings in the ex- 
ponent 1/3 in the general Boussinesq system conditionally, for slowly varying 
n. Another way by which the Nusselt number dependence on the Rayleigh 
number can be lower is if the Prandtl number is very high or in rotating (||14||) 
convection. If the Prandtl number is infinity then the upper bound is closer 
to 1/3. The equations of motion for infinite Prandtl number Rayleigh-Benard 
convection (1151) are 



Au + Vp = RTe z , V ■ u = 



(17) 



is time in- 



coupled with the advection-diffusion equation (g). Because (|T] 
dependent we say that T obeys an active scalar equation. 

An important observation, true even for the general case ([]]) is that, in 
view of the boundary conditions and incompressibility, not only the vertical 
component of velocity w but also its normal derivative ^ vanish at the 
vertical boundaries. Therefore we can write 

,/ \ T 7 ^ f Zl d 2 w . " , 7* ( d(T(-,a,t) -T(a,t))\ ' 

b(z,t) = ^ dzx j o —(., s ,t)ds j o [ g z ^-^Ua. (18) 



Consequently 



\b(z,t)\ < zi 



d 2 w(-,t) 

dz 2 



\\V(T-T)(;t)\\ 



(19) 



L°° 



loo is the sup-norm. One can express dr, in terms of T 



where 

Indeed, eliminating the pressure from ( |17D one has 

A 2 w = -RA h (T-T) 



T. 



(20) 



where Ah is the Laplacian in the horizontal directions x and y. Using the 
boundary conditions we may write this as 



w 



-RiAl^A^T-T) 



(21) 



where (A 2 DN ) 1 is the inverse bilaplacian with homogeneous Dirichlet and 
Neumann boundary conditions. Taking two z derivatives then yields 



d 2 w 
dz 2 



—RB(T - T) 



2 Wa 

Q z 2 \—DN) 



(22) 



(23) 



where the linear operator B is given by 

d 2 

The temperature equation obeys a maximum principle so that 

< T < 1 

holds pointwise in space and time. Therefore 

< |T_ T\ < 1 



holds. The operator B is not bounded in L°° but obeys a logarithmic extrap- 
olation estimate. This means that higher derivatives enter logarithmically in 
the bound; the estimate 



ll^(^ - T) ||^) < {l + log + i?} 4 



(24) 



follows from the bounds in (JTBl). The constant C\ can be computed ex- 
plicitly. Therefore, using the same kind of background as above in (|TTp we 
obtain 



N< / (T'(z)) 2 dz + 2 / £\t\z 



d 2 w(-,t) 
dz 2 



\\V(T-T)\\)dz (25) 
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and consequently 



N< — + 2E5l C X R [l + \og + R} 2 ^/E (26) 



5 

with C denned in (0), C\ coming from (H) and 



E= I 1 slp(s)ds. (27) 



optimizing in 5 < 1 we find 
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-i 



C 

and then letting P^swe deduce 



^^{l + log+i?} 2 ^ 



+ 1 



Theorem 2 There exists a constant C2 such that the Nusselt number for the 
infinite Prandtl number equation is bounded by 



4 

N < 1 + C 2 R^ {l + \og + R) 7 ns. 



If no additional information is given then, using n < N in the theorem 
above we recover the result 



N < l + C%R} {l + log + i?} 



of ( ||16|| ). The theorem brings in the exponent 2/7 as long as n is not varying 
too much with R. 

Discussion 

There is no proof that the Nusselt number cannot ever scale like R^ for ex- 
ceedingly large R. The experimental data do not seem however to encounter 
this behavior. Mathematically, the Nusselt number represents the maximum 
(among all possible invariant measures) of the expected value of the diameter 
of the global attractor in the energy dissipation norm. The functions on the 
attractor are not arbitrary, and may have certain properties that explain the 
experimentally observed bounds In this paper we showed that if the 

ratio 

2 



n (||V(T-r)|| 
N (\\VTf) 
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is small then the Nusselt number dependence on Rayleigh is depleted. 
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